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Transient Oscillation of Currents in Quantum Hall Effect of Bloch 
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We consider the quantum Hall effect of two-dimensional electrons with a periodic potential 
and study the time dependence of the Hall and longitudinal currents when the electric field 
is applied abruptly. We find that the currents oscillate in time with very large frequencies 
because of quantum fluctuation and the oscillations eventually vanish, for their amplitudes 
decay as \/t. 

KEYWORDS: integer quantum Hall effect, Chern number, linear response theory 



1. Introduction 

It is renowned that the Hall conductance in two dimensional electron systems under a 
strong magnetic field is quantized to an integer or a fraction multiplied by e^//i with very high 
accuracy.^ The relations between the conductance and topological numbers were discussed 
extensively,^"^ since the topological numbers take quantized values exactly. In the present 
paper, we discuss the integral quantization in noninteracting Bloch states. Thouless, Kohmoto, 
Nightingale and den Nijs (TKNN) showed using the Kubo formula that the quantized Hall 
conductivity is represented by the Chern number, which is a topological number defined on the 
two-dimensional torus (i.e., the magnetic Brillouin zone).^'^ The same result is also obtained 
from the adiabatic approximation. ^'^"^^ It would be an intriguing issue, at least from a purely 
theoretical point of view, that how the topologically quantized conductivity is modified when 
we go beyond the Kubo formula or the adiabatic approach. 

Interest in the TKNN theory was renewed recently in the field of ultra-cold atomic gases. 
The TKNN Hamiltonian is mapped to the Hamiltonian of a cold atomic gas trapped by a 
rotating optical lattice. Rotating Bose-Einstein condensates in a co-rotating optical lattice 
was indeed experimentally realized recently, which fueled the interest in the TKNN theory. 
The atomic gas system does not contain any perturbative effects coming from impurities or 
long range Coulomb type interactions. Hence, compared with the electron system in the solid 
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states, the atomic gas system is clean and the theoretical results of the TKNN theory can be 
applied without taking into account the corrections from such perturbations. An alternative 
method of applying an effective magnetic field to a cold atomic gas is also proposed. This 
method utilizes the internal degrees of freedom of cold atoms instead of the rotation of the 
system. The Hofstadter butterfly/^ which has been observed in a two-dimensional superlattice 
structure in a semiconductor heterojunction/''"^^ is predicted to be studied more easily using 
cold atomic gases. 

In this paper, we focus on the effect of a suddenly applied dc electric field on the integer 
quantum Hall effect of Bloch electrons.^'' The results are readily applied to the cold atomic 
gas trapped by a rotating optical lattice. We calculate the resulting current with the Kubo 
formula. ^^"^^ The linear response theory for an abruptly applied dc field was particularly 
investigated by Greenwood. We here follow Greenwood's formulation of the linear response 
theory. 

An interesting feature of our finding is an observation of fluctuation around the quantized 
conductivity, which is normally considered a very rigid quantity; we find that the Hall current 
has a time-dependent correction term to the Chern-number term in the TKNN theory. The 
Hall current jx and the longitudinal current jy oscillate in time with large frequencies because 
of quantum fluctuation, oscillation between different subbands. The oscillation eventually 
ceases and the time-dependent Hall current converges to the Chern-number term of the TKNN 
theory. The amplitude of the oscillation decays as 1/t. In the previous paper, ^'^ we already 
reported the existence of time-dependent correction terms. In the present paper, we present 
additional calculations particularly on the long-time behavior and on the time-dependent fields 
under an applied current. 

This paper is organized as follows. In §2, we derive the currents in the x and y directions 
following the Greenwood linear response theory. We derive the same results as in our previous 
paper, but under a different gauge. We also mention the correspondence between electron gases 
in a magnetic field and rotating cold atomic gases. In §3, we show that the time-dependent 
oscillation of the currents decays as 1/t and eventually ceases, and the Hall current approaches 
to a certain value obtained from the TKNN theory. Finally we give conclusions. In Appendix, 
we calculate electric fields under an applied current instead of currents under an applied field. 
We show that the voltages have similar time dependence. 

2. Time Dependence of Currents 

We consider noninteracting electrons in a periodic potential in the x-y plane. A magnetic 
field B is applied in the z direction. At time t = 0, we suddenly apply an electric field E{t) in 
the y direction. We calculate the currents of this system with the Kubo formula. The Kubo 
formula for a step-function external field is also known as the Greenwood linear response 
theory.^^ 
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Table I. Correspondence between Hamiltonians (1) and (4). 



Electron gas in a magnetic field 


Rotating cold atomic gas 






eB/2me 


n 


eEy 


Vy 



Using the Landau gauge, we write the Hamiltonian of the system as 

n = no-eyEyd{t), (1) 

where 

= ^[f'x + (P2/ + ei?x)VFp, (2) 

V, = f/ocos(^)+^osin(^). (3) 

In our previous paper, we treated the external field as a time-dependent vector potential. We 
here use the time-dependent scalar potential as in Eq. (1). We show below that the resulting 
formulae are the same. 

The Hamiltonian (1) can also describe a rotating dilute cold atomic gas trapped in an 
optical lattice.^^ To see this correspondence, let us consider a cold atom with mass ma confined 
in a harmonic potential. The periodic optical lattice which traps the cold atom rotates in the 
z direction with angular momentum fi. At t = 0, the optical lattice is tilted along the y-axis 
or accelerated in the y direction. The Hamiltonian of this system is written in the rotating 
frame as 

where = *(0, 0,^2) and L = p x x. We note that the centrifugal force is canceled because 
the frequency of the harmonic trap is the same as the frequency of the rotation, and the 
interactions between atoms are neglected. Equation (4) is also expressed as 

1 



n = -^ {p^ - mnyf + {py + mnx f +Vp-yVy9{t), (5) 
By substituting me, eB/2me, and eEy for ma, ^, and Vy, respectively (See Table I), we have 

Thus, moving from the Landau gauge to the symmetric gauge by the operator exp[i(eB /2h)yx], 
we see that the Hamiltonian (1) for an electron gas is identical to the Hamiltonian (4) for a 
cold atomic gas. 
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We consider the ratio (j) (= <&/$o) of the flux <I> (= Bab) per unit cell to the flux quanta 
$0 (= h/e). We put 

(7) 

where p and q are coprime integers. Because of the presence of the periodic potential, each 
Landau level splits into p sublevels. 

Let us first consider Hq. We write the eigenvalues and eigenfunctions of Hq as 

Wo \(t>Nm) = \<pNm) , (8) 

where the subscript N labels Landau levels and the subscript m labels sublevels in a Landau 
level {1 < m < p). We define the generalized crystal momentum hk in the magnetic Brillouin 
zone:^ < kx < 2TT/qa and Q < ky < 27r/6. Note that e'^^*" and e^^y^ are the eigenvalues of 
the translational operator. We define 

/Tq^ = e /To6 , 

\(l)Nm) = e'^-"" \uNm{k)) , (9) 

which satisfy 

^Ofe \uNm(.k)) = ENm{k) \uNm{k)) ■ (10) 

We thus block-diagonalized the Hamiltonian Hq into each subspace of k. 

Let us consider small Uq and treat the periodic potential as a perturbation in the subspace 
of a crystal momentum. Taking the lowest order terms into account, we obtain the wave 
function as^ 

UNmik;x,y) = ^ xn (x - qas - + kyi^ 

n=0 s=— oo ^ ^ 
X ^-ikx{x-qas-qan/p)^-2m(sp+n)y/b ^^-^-^ 

where £ = ^/h/eB is the cyclotron radius and Xn{x) satisfies 



dxXN{x) 



X 



2 2N + 1 



Xn{x). (12) 



£4 (2 

We note that UNm{k;x,y) in Eq. (11) satisfies the magnetic Bloch theorem: 

UNm{k; X + qa, y)e'^^^Py/'' = UNm{k; x, y) = UNm{k; x,y + h). (13) 

We have the eigenenergy within the perturbation as 

E^rnik) =hw^[N+^+ em{k), (14) 

where ujc is the cyclotron frequency. Here, e^(fe) and satisfy the following secular equation 
(the Harper equation) 



_wqb f 2Trq qbky\ ,„ Uo 
Uoe cos — -n - - — - dIL + -^e 
P P J 2 



^n+l^ikxqa/p _|_ ^n-l^-ikxqa/p 



em{k)d^. (15) 
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The coefficients satisfy (fm^ = and each Landau level splits into p subbands. 

We consider the currents caused by the electric field Ey9{t). In cold atomic gases, we 
can apply an effective electric field corresponding to Ey9[t) either by making use of the 
gravitational force tilting the harmonic potential^^ or by accelerating the optical lattice. We 
calculate the currents in the a (= x, y) direction in the form 

j^{t) = T,p{tf^, (16) 

where 



1 



-Px, 



me 



(17) 



Here p{t) is the density operator. 

Following Greenwood, we expand p{t) with respect to the electric field Ey and take the 
zeroth- and first-order terms into account: 

p{t)^po + piit). (18) 

The zeroth term po is the initial density operator, / Tre~^^° . With the help of the von 

Neumann equation for the density operator, pi is calculated as 

pi{t) = 4 r dt'ei^"(*'-*)/^'[-ey£;,^(0,Po]e-^^»(*'-*)/'^ 

l"' J— oo 



^ /* dt' ^-Ho^t'-m [-y, po] e-'^o^*'-*)/'^. 
Jo 



(19) 



We note that the lower bound of the integral on the second line of Eq. (19) is zero because 
of the step function in the perturbation, whereas the lower bound is negative infinity in the 
TKNN theory.^ By taking the trace in Eq. (16) with respect to the states in Eq. (11), we 
obtain the currents as 



jait) 



IVpi(t) 



E 



qab 

EE 



'^-^li'' dK f^^'^ dfc 



/f {ENrn{k)) 



l-Klqa Jo 27r/6 {ENm{k) - E^'m' {k)f 

X [{uNmm'^yk \uN'm'{k)) {uN'm'{k)\v^^ \uNm{k)) [l - ^-■<E^m{k)-E^,^,{k))t/n 
Nm N'm' -^MBZ 



X Im 



dko. 



'U'N'm' (k) 



duNmik) \ 

dka I 



-i{EKm{k)-E^,^,{k))t/h 
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where MBZ stands for the magnetic Brillouin zone. Here we used = e^'^-^He'^-^ and 



~\k x ik x ^ 



k 



(21) 



Furthermore, noting the relation Vy = [y,Tlo], we used 

{uNm{k)\ e-^^-*ye'^-* \uN'm'{k)) = Et,,^,{k)^- ENm{k) ^^^"^^^^l ""yk \^N'm'{k)) • (22) 

Let us put the Fermi energy in a finite gap between the moth and (mo + l)st subbands 
which belong to the lowest Landau level {N = 0). We consider the zero temperature. Hence 
the Fermi distribution satisfies fF{ENm) = 1 H N = and m < mo, and fpiENm) = 
otherwise. Thus we obtain 

-2Ey — - / — Im<^ 
27rn . MTK7. 27r 



Jait) 



+ 



E E 

m<mo m'>mo 



MBZ 



dk. 



y 



UOm' (k) ) ( Uom' (fe) 



dkn 



1 _ Q-'^{Eom{k)-EQ^,{k))t/h 
duQ„i{k 



E E 

m<mo N'>l,m 



dk,. 



(fe) 



Oka 



1 _ Q-i{Eom(.k)-Ej^,^,ik))t/h 



-2Ey-— / 7i-Im<! 
27rn 7mbz 27r 



m<.'rrn) ■m'>'m{) 



\ 9ky 



UOm'{k) ){ UQm'{k) 



dkn 



i{Eo^{k)-Ea^,{k))t/h 



1 



m<mo 

where we used the fact that for iV' > 1, we have 



i{Eom{k)-Eim{k))t/h 



duQmjk) 
dkn 



(k) 



duomik) 



(23) 



(24) 



Oka I 2 

We can see that the time dependence of the current is due to the quantum fluctuations or 
quantum oscillations between various sets of discrete levels. We ignore the quantum fluctuation 
between ^o^(fc) and Eijn{k) because its frequency, which is proportional to Eojn{k) — Eijn{k) , 
is very large compared to the frequency of the fluctuation between different subbands of the 
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lowest Landau level, which is proportional to EQjn{k) — E^^iik). Thus we obtain 

m) = ^[Nch + Aa,{t)], (25) 

jyit) = (26) 



where 



Ch 



m<mo 



MBZ 27ri 



dk-r 



dk„ 



c.c. 



(27) 



Aay{t) 



m<mo m'>mo 

\ dky 

Y. Y. 

m<mo m'>mo 



MBZ TI" 
fk 

MBZ 



/ duomjk) 
\ dky 



dk^ 



Uom'{k) 



Q^{^n^>ik)-em{k))t/h 



(28) 



sin[(e„/(/c) - em{k))t/h] . 



(29) 



Note that A'^ch is the Chern number and takes integer values.^' ^ The time-dependent correction 
terms A.ax{t) and Aay(t) express quantum fluctuation between different subbands of the 
lowest Landau level. These arc expressed as the sum of different oscillating modes whose 
frequencies are determined by the energy difference €m'ik) — e„j(fc). 

Hereafter, we show results of numerical calculation of the currents jxit) and jy{t). Nu- 
merical calculation is carried out in a way similar to the Kubo formula for a dc field:^^'^^ the 
integrals in Eqs. (27), (28), and (29) are performed with random sampling of kx and ky. In 
the calculation, we set a = b. We here consider, for example, the following three cases: (i) 
p/q = 5/4 and tjiq = 2 (A^ch = 2), (ii) p/q = 7/6 and mo = 3 (A^ch = 3), and (iii) p/q = 7/6 
and mo = 1 (A'^ch = !)• 

The band structure in the case (i) is shown in Fig. 1. In the figure, the Fermi energy 
that we choose is plotted with the dashed line. Figure 2 shows the currents jx{t)h/e'^ Ey and 
jy{t)h/e'^Ey in the case (i). In the calculation, we put Uq = O.lmeV and a = b = lOOnm 
as tipical values for quantum Hall systems on a semiconductor heterojunction. The currents 
oscillate irregularly reflecting the fact that the energy spectra em=2ik) and em=3{k) in Fig. 1 
strongly depend on k, and Aax{t) and Aay{t) are written as the sum of sinusoidal functions 
with different frequencies (Eqs. (28) and (29)). The insets show the long-time behavior of the 
currents. As we show in the next section, Ao"^(t) and A(Ty{t) vanish for large t. 

Similarly, the band structure in the cases (ii) and (iii) is shown in Fig. 3. Since Nq\i in 
Eq. (27) depends on mo, A'ch changes when we change the Fermi energy. In the figure, the 
Fermi energy for the case (ii) is plotted with the dashed line and that for the case (iii) is plotted 
with the dotted line. Figure 4 shows the currents jx{t)h/e^Ey and jy{t)h/e^ Ey in the case 
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27t/qa 




Fig. 1. Case (i): The band structure of as functions of kx (left) and ky (right). The flux ratio 
p/q = 5/4. We place the Fermi energy (the dashed line) between the second and third subbands. 




0.5 1 1.5 2 0.5 1 1.5 2 

t [ns] t [ns] 

Fig. 2. Case (i): The currents jx{t)h/ Ey and iy{t)h/ e^Ey arc shown as functions of time. We set 
p/q = 5/4, mo = 2, Uq = O.lmcV, and a — b — lOOnm. The dashed lines show the convergent 
values of jx{t)h/ e^Ey and jy{t)h/e^Ey {Nch (= 2) and 0, respectively). The insets show long-time 
behaviors of jx{t)h/e^Ey and jy{t)h/e^Ey. 



(ii). The CTirrcnts oscillate irregularly because of contributions from different frequencies. The 
insets show the long-time behavior of the currents. Figure 5 shows the currents jx{t)h/e'^Ey 
and jy{t)h/ Ey in the case (iii). In this case, the currents oscillate rather regularly because 
the first and second bands in Fig. 3 are almost flat, and /S.(Tx{t) and Acry(t) are almost 
monochromatic. The insets show the long-time behavior of the currents. 

We remark the following three points. Firstly, the currents jx{t) and jy{t) are gauge 
invariant. We can also obtain the same results by using the time-dependent vector potential as 
we did in the previous paper. Secondly, if a dc current instead of a voltage is abruptly turned 
on, the voltages in the x and y directions temporarily vary in the same manner as Eqs. (25) and 
(26), i.e., the period of the oscillation is given by the energy difference between two sublevels 
(see Appendix). Finally, although the electric field is given by the step function here, we axe 
able to calculate the time dependence of the currents for an arbitrarily time-dependent electric 
field by following the machinery of the Kubo formula (See Eq. (A-2) below). 
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Fig. 3. Cases (ii) and (iii): The band structure of as functions of (left) and ky (right). The flux 
ratio p/q = 7/6. We place the Fermi energy for the case (ii) between the third and forth subbands 
(the dashed line), and the Fermi energy for the case (iii) between the first and second subbands 
(the dotted line). 
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Fig. 4. Case (ii): The currents jx{t)h/e^Ey and jy{t)h/ e^Ey are shown as functions of time. We set 
p/q = 7/6, mo = 3, ?7o = O.lmeV, and a = b = lOOnm. The dashed lines show the convergent 
values of jx{t)h/e^Ey and jy{t)h/e'^Ey {Nch (= 3) and 0, respectively). The insets show long-time 
behaviors of jx{t)h/e^Ey and jy{t)h/e^Ey. 



3. Long-Time Behavior of Currents 

Let us study the currents after a long time. We show that |A(T^(t)| and |A(7jy(t)| decay 
as 1/t using the Riemann-Lebesgue theorem:^° limt^oo J^'' g{uj)e^^duj = 0, where ^(a;) is 
uniformly convergent. 

Both Aax{t) and Aay{t) are expressed as {a = x, y) 

Aaa{t)= [ d^klmgl^l,{k)e<'rn'ik)-.m{k))t/!i_ ^^^^ 

m<mo m'>mo "^^^^Z 

We define 

— ^m'jk) — e.m{k) , , 

'^mm'(«j = ^ , (^-Ij 

f^O = <^mm'{kx = 0, ky), LV^iqa = ^mm'{kx = T^lm, ky). (32) 
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Fig. 5. Case (iii): The currents jx{t)h/e'^Ey and jy{t)h/e'^Ey are shown as functions of time. We set 
p/q = 7/6, mo = 1, Uq = O.lmeV, and a ^ b = lOOnm. The dashed lines show the convergent 
values of jx{t)h/e^Ey and jy{t)h/e^Ey {Nch (= 1) and 0, respectively). The insets show long-time 
behaviors of jx{t)h/e^Ey and jy{t)h/e^Ey. 



dkr 



Hence, 

Aaa{t) = 2 ^ ^ / dkyl du, 

m<mo m'>mo ^ 
m<mo m'>mo •''^^ 

where uja = min (wq, a;^/ga) and iOh = max (wq, a;^/qa)- 

We note that J"^'' dwe"^* = (e^'^''* — e''^"*)/it and therefore this integral decays as Ijt. We 
express 5^^, (a;) as 



(33) 



(34) 



where u> and t;> are positive and n< and f< arc negative m. LOa < < ^b- We write the 
maximum and minimum of these functions as u™^^ = max|u>(u;)|, u™"^ = min|u>(a;)|, etc. 
Then we see the integrals, da;n"""e''^*, duju^^^c^^^ , etc. also decay as 1/t. We have 



< 



< 



/ ' 



da;^TO''u>(a;)e" 



dUmm'U<{ljj)e^'- 



< 



< 



/ ' 



du iw^^/t 



etc. Therefore 



Thus we have shown 



/ ' 



du;mm'£?^i'(w)e"^"""' 



(35) 
(36) 

(37) 



|A(Ta(i)| ~ {a = x,y). 



(38) 



In Figs. 6, 7, and 8, we show logarithmic plots of |A(Ta;(t)| and |Acr2^(f)| in the three cases 
(i) p/q = 5/4 and mo = 2 {Nch = 2), (ii) p/q = 7/6 and mo = 3 (A^ch = 3), and (iii) 
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0.001 




Fig. 6. Case (i): Logarithmic plots of the long-time behavior of |A(T2;(t)| and |A(Tj,(t)|. We also draw 
the dashed lines 0.4/< (left) and 0.6/t (right) to see |Acra;(t)| and \Aay{t)\ decay as 1/t. The 
parameters are the same as in Fig. 2. 



p/q = 7/6 and rriQ = 1 (Nch = !)• In all cases, \Aax{t)\ and \Aay{t)\ indeed decay as 1/t. 
Thus, the response of the system to the temporal change of the external field disappears in 
nano-second order even if there is no dissipative mechanism. 

Since the correction terms Aax{t) and Aay{t) decay as 1/t, in the limit t — > oo, we obtain 

^^^'^^^ = iS'^C^' 

jy{t ^oc) = 0. (39) 

This Hall current was first obtained by Thoulcss et al? 

When the bands em(fc) and em'{k) are nearly flat as in the case (ill), we can explicitly 
calculate the time dependence of Aaa{t). In this case, ui, — uia is very small and Eq. (33) can 
be approximated as 

Aa^{t)^Y. El°^/ ^ky-gl^liky) da;„^,e-wt, (40) 

m<mo m'>mo ''^ •' 

where g^li,{ky) 9^^' i^^mm' , kyj^ is independent of Umm'- We note that 

P da;e-* = ^ sin e^K+-.)V2. (41) 

This integral decays as 1/t and its amplitude has two kinds of oscillations. The period of 
one oscillation is inversely proportional to — Ua and the period of the other is inversely 
proportional to LOb + coa- We note that the difference LOb — coa is very small and cOa — cob- 
Therefore, the frequency of exp[i{L0b + u!a)t/2] is given by the energy difference between e^(fe) 
and e^/(fc), and the period of the beat 4Tr/{ub — coa) is very long. The 1/t decay is revealed 
for a time longer than the period of the beat. 

Thus, in nearly fiat-band cases, it is easier to observe the 1/t dependence because |A(7j;(t)| 
and |Acrj^(t)| decay rather slowly and survive for a long time as is seen in Fig. 8. 
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Fig. 7. Case (ii): Logarithmic plots of the long-time behavior of \Aax{t)\ and |Ao'j,(t)|. We also draw 
the dashed lines 1.0/t (left) and 1.5/t (right) to see \Aax{t)\ and \Aay{t)\ decay as 1/t. The 
parameters are the same as in Fig. 4. 
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0.0001 



Fig. 8. Case (ii): Logarithmic plots of the long-time behavior of \Aax{t)\ and \A(Ty{t)\. We also draw 
the dashed lines 2.5/t (left) and 2.5/t (right) to see \Aax{t)\ and \Aay{t)\ decay as 1/t. The 
parameters are the same as in Fig. 5. 



4. Conclusions 

Using the Greenwood linear response theory, we studied the time dependence of the cur- 
rents in the quantum Hall effect when the electric field is suddenly turned on. We found that 
both jx (t) and jy {t) oscillate because of the quantum fluctuation between two subbands which 
straddle the Fermi energy. These oscillations decay as 1/i and eventually cease. In the limit 
t —>■ oo, jx{t — > oo) is given as the Chern number A^ch multiplied by e^/h as Thouless et al? 
obtained. As is discussed in Appendix, the electric flelds oscillate in time in the same way as 
jx{t) and jy{t) when, in reverse, the current is applied abruptly at t = 0. 

We showed that the ratio of the Hall current and the suddenly applied dc field is decom- 
posed into the sum of a constant term and a time-dependent term (Eq. (25)). The constant 
term is the conductivity for the dc field applied for infinite time and given by the Chern 
number Nq^. Thus, the time-dependent term /^.a^it) expresses a correction to the Chern 
number term. In other words, ^CTxit) can be regarded as the fiuctuation around the Chern 
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number. It will be remarkable to observe the fluctuation experimentally, since the quantiza- 
tion to the Chern number is normally regarded as very rigid. This fluctuation, which stems 
from transitions between different subbands, decays as 

Thus, the response of the system to the temporal change of the external field decays as 1/t 
even if there is no dissipative mechanism. The amplitude of the decay gets large if the bands 
that give large contribution to A(7^(t) are nearly flat. In this case, the 1/t-decay survives for a 
long time. In a quantum Hall system on a semiconductor heterojunction, this power-law decay 
of the order of nano-second might be difficult to observe experimentally because relaxation 
time due to impurity scattering, etc. is of pico-sccond order. Cold atomic systems under an 
artificial magnetic field may overcome these difficulties. In experiments of a Rubidium cold 
atomic gas trapped by a rotating optical lattice, the time scale of the power-law decay is of 
the order of millisecond for a = 5 ~ l/nm, Uq ~ O.lneV, and large O so that p/q ^ 1. 
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Appendix: Measuring Electric Fields under an Applied Current 

Here, we calculate the voltage for the applied dc current that is switched on abruptly. This 
situation matches current-controlled experiments. (For theoretical reasons, in the main body 
of the paper, we calculate the current under the applied voltage.) We find that the voltage also 
oscillates. Both temporal oscillations of the current and voltage are caused by the quantum 
fluctuation between two subbands. 

We apply the current suddenly at t = in the y direction, jxit) = 0, jy{t) = Jy9{t), and 
obtain Ex{t) and Ey{t). This may be closer to the experimental situation. 

Since we assume linear response, we have (a = re, y) 





By Fourier transform, we obtain 




(A-2) 



I3=x,y 
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where Jo ("J^) = / Ja(i)e^'^*dt, etc. We put an infinitesimally small r/ > to ensure the causality: 
o'aisit) = for i < 0. We define the resistivity as 

/3=x,y 

where 

-(^xy{^' + ill) 



Pyyiuj + ir]) 



Therefore we obtain electric fields as 



a^y{u + \ri) 



(A-4) 

(A-5) 
(A-6) 

(A-7) 

uj + \ri ^ ' 

Wc obtain the conductivities dxy{ijo) and dyy{u)) with the help of the calculation in §2. We 

first note that 



1 

2vr J ""^ aly{uj + ir/) + aly{u: + i??) 

CFyy{io + i??) 



~jyi^>-'''\ 



where 



Ey{u) = j AtEye{t)S'^+''''^^ 

Using Eqs. (25) and (26), we have 



iE, 



U} + 17] 



axy{uj + ir]) 



Ey{u) 

„2 



27rin 



(a; + irf) 



OJ + 17] 



X 



MBZ 27r \ dky 
1 

u + \r] + {em'{k) - em{k)) /h 



+ ( E E - E E 

^m<mo m'>mo m>mo m'<mo 



dkor 



ayy{uj + ir]) 



UJ 



Ey{u) 



(A-8) 



(A-9) 



X 



27rin 



m<mo m'>mo 'rn>mo m' <mo j 



d^fc / duo„,,{k) 
MBZ 27r \ (9A;. 



2/ 



LO + ir] + (e^'(fc) - em(fc)) 



(A-IO) 
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By plugging Eqs. (A-9) and (A- 10) into Eqs. (A-5) and (A-6), we obtain the electric fields: 



2'Kh 



s 



m<mo m'>nn) 



duom{k) \i(,^,(k)-em(k))t/h^ (A-ii) 

dkx / 



Ey{t) = 



2Trh 



c2 '^y 



d^fc g ( era'{k) - era{k) 

m<mo m'>mo 



h 



\ 9ky 



UOm' {k) 



sin[(e^/(fc) - emik))t/h)] , 



where 



2'Kh 



(A-12) 



(A-13) 



Note that Ex{t) and Ey{t) have the same time dependence as jx{t) and jy{t) in §2; the period 
of the oscillation is dominantly given by the energy difference between two subbands which 
straddle the Fermi energy. 
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